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CP nonconservation in the leptonic sector 



Petre Di^Cl 

In this paper we use an exact method to impose unitarity on moduli of neutrino PMNS matrix 
recently determined, and show how one could obtain information on CP nonconservation from a 
limited experimental information. One suggests a novel type of global fit by expressing all theoretical 
quantities in terms of convention independent parameters: the Jarlskog invariant J and the moduli 
\U a i\, able to resolve the positivity problem of \U e 3\- In this way the fit will directly provide a value 
for J, and if it is different from zero it will prove the existence of CP violation in the available 
experimental data. If the best fit result, \U e 3\ 2 < 0, from M. Maltoni et al, [New J.Phys. 6 (2004) 
122] is confirmed, it will imply a new physics in the leptonic sector. 
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PACS numbers: 12.15.Hh, 12.15.-y, 14.60.Pq 



I. INTRODUCTION 

Neutrino experiments have shown that neutrinos have 
mass and oscillate, the last property suggesting a mixing 
of leptons similar to that of quarks for more ref- 

erences on experimental data see, e.g" QL'Q. The pres- 
ence of the lepton mixing opens the possibility to see a 
CP violation in the leptonic sector, phenomenon which is 
considered an indispensable ingredient that could gener- 
ate the excess of the baryon number of Universe 0] , since 
the baryon number can be generated by leptogenesis Q . 
Excluding the LSND anomaly, see [9l]-[l0j|. all neutrino 
data are explained by three flavor neutrino oscillations, 
and the determination of oscillation parameters is one of 
the main goals of phenomenological analyzes. The LSDN 
signal for —¥ T> e has been interpreted as a possible exs- 
tence of a fourth light neutrino, or as a signal of CPT 
violation. The current opinion is that neither explana- 
tion provides a good fit to all existing neutrino data, see 
e -g- UM- Recently two detailed combined analyzes of all 
neutrino data have been published, see Refs. [5-6], that 
materialized in upper and lower bounds on sin fly, these 
bounds being converted in Ref. [12j j into intervals for the 
moduli of the PMNS unitary matrix [l3| . In both above 
papers the CP phase 6 was not considered a free param- 
eter, e.g., only the particular cases cos<5 = ±1 have been 
considered in Q. We shall write the numerical moduli 
matrix as 



V = 



0.835 ±0.045 
0.355 ±0.165 
0.365 ±0.165 



0.54 ±0.07 
0.575 ±0.155 
0.590 ±0.15 



< 0.2 
0.7 ±0.12 
0.685 ±0.125 



(1) 



where the central values are given by half the sums of 
lower and upper bounds entering [l2| , and the errors cor- 
respond to a three sigma level, excepting V e s for which 
only an upper bound is given. Complementary infor- 
mation concerns the absolute value of the two neutrino 
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mass-squared differences [j| , but we have no information 



• the magnitude of £/ e 3-element of the leptonic mixing 
matrix, or the value of the mixing angle #13; 

• the existence of leptonic CP violation, or a value for 
the Dirac phase S 7^ 0, tt entering the PMNS matrix [l3j ; 

• the order of the mass spectrum, or the sign of Am^; 

• the magnitude of neutrino masses, etc.; see, e.g., [bH ]- 

Although it is believed that only the next generation 
experiments will give an answer to the first above items, 
[l7j, we consider that by using unitarity we can obtain 
some information from experimental data as given in ([T]) 
for both the existence of CP violation in leptonic sector, 
and the parameter #13, all the more as measurement of 
#13 is considered a first mandatory ingredient for inves- 
tigation of CP leptonic violation in — > v e transitions 
and for the mass hierarchy determination. Concerning 
CP violating phase, 5, the general opinion is that it will 
require a major experimental effort because of the intrin- 
sic difficulty to measure it, see e.g. Ref. fTs| ] - 

The current phenomenological analyzes of data are 
given in terms of mixing angles 9ij, see, e.g., |5j-[6|, and 
phase S which all depend on the chosen parametriza- 
tion of the unitary matrix one starts with. In the pa- 
per we propose another method for analyzing the data 
by expressing all measurable quantities directly related 
to oscillation amplitudes in terms of moduli, |f/ a j|, and 
Jarlskog invariant J, [l9j . From a theoretical point of 
view this approach has the advantage that the recon- 
struction of a unitary matrix from knowledge of moduli 
leads essentially to a unique solution, see, e.g., Ref. [2fj| . 
From an experimental point of view the method allows 
to get from nowadays limited information on neutrino 
physics all moduli of the unitary matrix such that one 
can predict the approximate form for the oscillation am- 
plitudes, for example v T — > which could be of great 
help for experimenters. It is easily seen that the experi- 
menters have measured quantities directly related to os- 
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dilation amplitudes, which contain information on the 
moduli \U e i\, \U e2 \, | C/^s | , and \U t3 \. The above mod- 
uli are correlated by unitarity to | U e3 \ 

\u e3 \ 2 = i - \u el \ 2 - \u e2 \ 2 = i - \u^\ 2 - \u T3 \ 2 (2) 

By using ((2|), from the first row and from the third col- 
umn in (HJ) one gets the bounds 

-0.1465 < |C/ e3 | 2 < 0.155, 
-0.3285 < \U e3 \ 2 < 0.35 

The moduli of lower and upper bounds in the above 
relations are of the same order of magnitude, which 
means that the possibility V 2 3 < is not excluded. In 
fact the best fit of the KamLAND data in Ref. [5[ gave 
V 2 3 = —0.13; see also Figure A. 2 in Ref. @ and, respec- 
tively Figure 14 in Ref. 6]. Hence the phenomenologists 
have to take into account this possibility, which, if con- 
firmed, could be interpreted as a breaking of the estab- 
lished rules, signaling really a new physics. A slight mod- 
ification of the theory can take it properly into account. 
After nonvanishing masses, such a surprise coming from 
neutrino physics is not a priori excluded. 

As we will show in the paper our approach makes full 
use of unitarity property, which implies many correla- 
tions between all moduli of a unitary matrix. For a 
complete determination of an (exact) unitary matrix the 
knowledge of four independent moduli is necessary, see, 



e.g., Ref. [20|. Taking into account that a 3 x 3 unitary 
matrix has nine moduli, the choice of four ones is com- 
pletely at our disposal. Hence for analyzes we may use 
moduli which, apparently, do not appear in the measured 
oscillation amplitudes. In fact there are 57 groups of in- 
dependent moduli, and this fact, coupled with unitarity 
constraints, gives us an opportunity to find all moduli, 
even nowadays we have only a limited experimental in- 
formation. 

In the present paper we explore the possibility of ex- 
tracting information on CP violation in the leptonic sec- 
tor starting with the current knowledge on moduli matrix 
([1]). We find that within the moduli intervals given by 
((T|) there are values corridors that are compatible to the 
existence of at least one unitary matrix, in fact one can 
construct a continuum of solutions, a result which could 
be interpreted as a signal of CP violation in the leptonic 
sector. A clear signal of CP violation will be a value 
,/ 7^ for the Jarlskog invariant, a value which can be 
provided by a global fit suggested in the paper. If the 
new fit will confirm the result, |J7 e 3| 2 < 0, by Maltoni et 
al, this will be equivalent with the existence of new 
physics in the leptonic sector. 

The theoretical tool that will be used in the following 
is unitarity of the PMNS matrix, and the constraints 
on physical observables generated by this property. We 
make full use of phase invariance property to write the 
matrix under the form 



C12C13 C13S12 S13 



U = -c 23 s 12 e ld - c 12 s 23 s 13 c 12 c 23 e id - s 12 s 23 s 13 s 23 c 13 (4) 



si 2 s 23 e lS - ci 2 c 23 s 13 -ci 2 s 23 e lS - si 2 c 23 si 3 c 23 ci 3 



and the notation is Cij — cos 9ij and 
9ij are the mixing angles with ij = 12, 13, 23, and <5 is 
the Dirac phase which encodes CP violation. 

The matrix U is assumed and built as a unitary matrix, 
but nobody assures that the moduli | Uij | , the angles 9ij , 
or the phase 8 extracted from experiments come from a 
unitary matrix. Usually it is assumed that irrespective 
of how the measured data are, they are compatible to 
the existence of a unitary matrix. A problem which was 
not considered in physical literature is the property of 
unitary matrices to be naturally embedded in a larger 
class of matrices, the set of double stochastic matrices. 

A 3 x 3 matrix M is said to be double stochastic if its 



sin^y, where are double stochastic, matrices from the set 



entries satisfy the relations [21 1 



> 0, Ei=i ™>ij = 1, 



(5) 



All such matrices form a convex set, i.e. if Mi and M 2 



M 3 = x Mi + (1 - x) M 2 , with < x < 1 (6) 

are double stochastic [2l[. The unitary matrices are a 
subset of this larger set if we define the entries rriij of the 
double stochastic matrix by relations 



Mj — \Uij\ 



(7) 



One easily sees that relations ([S]) are satisfied because of 
the unitarity property 



(8) 



where f denotes the hermitian conjugate matrix, and -Z3 
is the 3-dimensional unit matrix. The subset of dou- 
ble stochastic matrices coming from unitary matrices is 
known as unistochastic matrices. It is well known in the 
mathematical literature that there are double stochastic 
matrices that do not come from unitary matrices, see, 
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e -g-; [HI - 01- This unpleasant feature rises a few prob- 
lems which have to be solved before trying to do a fit of 
the experimental data. The first one is to find the neces- 
sary and sufficient conditions the data have to satisfy in 
order they should come from a unitary matrix. In other 
words we have to find the separation criteria between the 
two sets. After that we have to provide an algorithm for 
the reconstruction of U from the given data. With these 
aims in view we have developed a formalism for express- 
ing the necessary and sufficient conditions for a 3 x 3 
double stochastic matrix to come from a unitary one. 
On the other hand the embedding relation ([7]) suggests 
that a "natural" parametrization of 3 x 3 unitary matrix 
could be that given by four independent moduli. This can 
be done with a proviso: cos 5 as the function of moduli 
must take physical values, i.e. — 1 < cos<5 < 1. The use 
of moduli is also justified by the fact that mixing angles 
and CP phase are not well defined objects, they depend 
on the form of the unitary matrix. In this respect there 
are at least two popular forms: the PDG proposal (25$ . 
and the classical proposal (26$ by Kobayashi-Maskawa, 
and the physics has to be invariant, i.e. not dependent 
of one or another choice. The moduli are such invari- 
ants. The above statements can be formalized by the 
following (obvious) axiom: irrespective of the physical 
processes where the moduli appear, or the phenomeno- 
logical models used in analyzing the data, the numerical 
moduli values should be the same. 

We used the novel formalism to develop an algorithm 
for the reconstruction of a unitary matrix from a double 
stochastic one when it is compatible to a unistochastic 
one. In this sense we modified the standard x 2_ test to 
allow the implementation of all unitarity constraints, for 
a better processing of experimental error affected data. 

To see that the found unitarity constraints are nec- 
essary to be used in any phenomenological analysis we 
constructed an exact double stochastic matrix, from the 
moduli ([1]) , which is perfectly acceptable from an experi- 
mental point of view, i.e. all the moduli are very close to 
the central values in (TTJ) being within the experimental 
error bars, but which is not compatible with a unitary 
matrix. 

Last but not least we used the convexity property of 
double stochastic matrices to develop a method for do- 
ing statistics on unitary matrices through their squared 
moduli, problem which was an open one in the literature 
until now, see, e.g., Refs. [13]- [H[ ■ 

The paper has the following structure. In Sec. II we 
present our phenomenological model that makes a novel 
use of the unitarity property of the PMNS matrix, the 
main point being the focus on the condition that sepa- 
rates double stochastic matrices from unitary ones. This 
condition is, in the same time, the consistency condition 
between the theoretical model, represented by the PMNS 
matrix, and the experimental data, as those appearing in 
(fTJ). In Sec. Ill we present the reconstruction algorithm 
of unitary matrices from experimental data, when they 
are compatible, and show how errors affected data can 



be consistently included in the formalism. We propose a 
two terms x 2 -test which has to contain a piece enforcing 
the fulfillment of unitarity constraints, and another piece 
which will properly take into account experimental data. 
In Sec. IV we propose a new type of global fit, done in 
terms of invariant parameters, J and moduli |fy|, which 
can resolve the existence of CP nonconservation in the 
leptonic sector by using the current experimental data, 
or, alternatively, prove the existence of new physics in the 
leptonic sector. In particular we show that around central 
values from data ([T]) there is a continuous (approximate) 
unitary set of matrices, all of them being consistent with 
S fa 7r/2, for a broad range of numerical V e3 values. By 
changing the moduli values around one looks for unitary 
matrices one gets other values for S and V e3 . The paper 
ends by conclusions. 



II. PHENOMENOLOGICAL MODEL 

By phenomenological model we will understand in the 
following a relationship between the entries of a double 
stochastic matrix, provided in general by experiments, 
and the entries of a unitary matrix. The experimental 
data will be the numbers affected by errors entering fl}. 
We remind that in order to find allowed ranges for the 
moduli \Uij\, as those given in (TT]), the solar neutrino and 
atmospheric neutrino measurements, as well as the reac- 
tor neutrino experiments have been used, see (EL 1^. fl2l] . 
The analyses lead to numbers for all moduli, but \U e3 \, 
from the first row and the last column. Unfortunately 
only three moduli are independent, because the unitar- 
ity gives 

\U el \ 2 + \U e2 \ 2 = \U^ 3 \ 2 + \U T3 \ 2 (9) 

Taking into account the embedding ([7]), which does not 
imply the unitarity fulfillment by experimental data, we 
will use different notation for the theoretical quantities, 
and for experimental data. Hence in the following we as- 
sume that the data are contained in a numerical matrix, 
V, whose entries are positive 





( V e i 


v e2 


v e3 


V = 










{ Vrl 


V T2 


v r3 



and the current experimental data on the moduli V a t , a = 
e, n, t, i = 1, 2, 3, are summarized in matrix ([l}. From a 
phenomenological point of view the main problem to be 
solved is to see if from a numerical matrix as (fTJ) one can 
reconstruct a unitary matrix as ([3]). For that we define a 
phenomenological model as follows 

V 2 = \U\ 2 (11) 
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relation that has to be understood as working entrywise 
leading to the following equations 

V-2 2 2 T/2 2 2 T/2 2 

el — c 12 c 13> V e2 — s 12 c 13, v e3 ~ S 13 

V-2 2 2 T/"-^ ^ ^ 

Ai3 — s 23 c 13> "tZ — c 13 c 23> 

V'tl = Sl2c| 3 + «13S23 C 12 + 2S12S13S23C12C23 COS 8, 

V^ 2 = C12C23 + Si 2 S^ 3 S2 3 - 2si2Si 3 S23Cl 2 C23COs5, (12) 

2 2 2 2 2 2 

V T l = S13C12C23 + S12S23 - 2S12S13S23C12C23 COS(5, 

2 2 2 2 2 2 

V t2 = S12S13C23 + C12S23 + 2S12S13S23C12C23 COS (5 

together with the double stochasticity relations 

i=l,2,3 i=l,2,3 

E l^i| 2 -l= E V&-1 = 0, J = 1,2,3 (13) 

The above phenomenological model is similar to that 
proposed by Wolfenstein, [29|, i.e. it is a direct relation- 
ship between the measured values Vij and the theoreti- 
cal parameters Sij and S entering U. The difference be- 
tween the two approaches is the explicit use of the double 
stochasticity properties (|13j) . and making no approxima- 
tions on right hand side Eqs. (HU). Also important is 
the fact that on left hand side of Eqs. (|12p we have sets 
of nine numbers, Vij, which all could be obtained from 
experiments, and on right hand side there are four inde- 
pendent parameters. Hence checking the consistency of 
the system of Eqs. (IT2t is a natural problem and it has 
to be resolved. Since the matrix U is parameterized by 
four independent parameters, e.g. Sjj and S, it is not at 
all obvious that the nine equations (fT2")l have a physical 
solution for all experimental values allowed in ([I]). On 
the other hand we have to be aware that the parameters 
and S are convention dependent and have no intrinsic 
theoretical significance, and this fact rise the problem of 
using only invariant functions generated by the entries 
of the matrix (j4j. Such functions are the Jarlskog in- 
variant J and the moduli V a i, see e.g. [30l |. and this is 
the main reason for starting with moduli as independent 
parameters in our approach. 

We said before that for n > 3 there are double stochas- 
tic matrices which are not unistochastic, see, e.g. [2l| - 
[24| , hence a novel problem arise: what are the neces- 
sary and sufficient conditions for their separation? In 
this resp ect there are only two mathematical references, 
[23|-|23|. where theorems are given on the necessary and 
sufficient conditions for a double stochastic matrix to be 
also unistochastic. See also Refs. [3lj]-[32| for the physi- 
cists point of view. The theorems are only existence the- 
orems, i.e. they do not provide a constructive method 
for recovering a unitary matrix from a double stochastic 
one when they are compatible. However they are impor- 
tant since show the equivalence between unitarity and 
the positivity of (any) area generated by the so called 
unitarity triangles, which is one form of the consistency 
condition between double stochastic matrices and uni- 
tary ones. Here, by using the phenomenological model, 
Eqs. pT1) - (|13l) . we provide an alternative method for the 



separation criteria between the two sets, the consequence 
being a new form of the consistency conditions between 
the experimental data and the theoretical model, which 
leads eventually to a constructive algorithm for recovery 
of unitary matrices from error affected data. 

To see how the phenomenological model defined by the 
relations ([TT]) -(fl"3" |) works, let simplify the things and as- 
sume for a moment that the relations (1131) are exactly 
satisfied by the experimental data, implying the data are 
not affected by errors. Then it is an easy matter to find 
from the first five relations (fl~2|) three independent ones 
which give a unique solution for the Sij, ij — 12, 13,23. 
In other words, if the experimental numbers satisfy the 
relations 

K 2 i + K 2 2 + v e % = 1 , v e % + v 2 3 + v 2 3 = i (14) 

we get always a physical solution for which is unique 
under the condition < < 1, and depends on the 
three chosen independent moduli. Substituting this so- 
lution in the last equations one gets four equations for 
cos 5, that lead to a unique solution for it. But nobody 
guarantees that the solution will satisfy the physical con- 
straint 

-1<cos<5<1 (15) 

Thus the necessary and sufficient conditions that guaran- 
tee the consistency between the data and the theoretical 
model are given by: 

0<Sij<l, ij = 12, 13, 23, -l<cos<S<l (16) 

The constraints (|T6]) together with relations (|12[) - (fT3)) 
prove the unitarity property of data, and provide an al- 
gorithm for recovery of unitary matrices from their ex- 
perimental moduli. 

The embedding relations Q of unistochastic matrices 
into the double stochastic set show that the number of 
independent parameters is the same for the two sets, i.e. 
in our case it is equal to four. Now we construct a double 
stochastic matrix by using the constraints (fT3|) and the 
entries from the data set, (f]~0|) . In order to simplify the 
formulas form we make the notation 

V e i = a, V e 2 = b, Vfti = c, Vfja = d, 

V^ = e, V Tl =g, V T2 = h, V t3 = f (17) 

As we already remarked the measured parameters a, b, e, 
and / are not independent since unitarity implies the 
relation a 2 + b 2 = e 2 + f 2 . Thus in principle we have 
to chose only three of them. Our choice for the first two 
matrices is, a and 6, since they have the smallest errors, 
and c and e, respectively, d and /. The third matrix will 
be generated by c, d, g, h, i.e. by moduli that are not yet 
measured. If the independent parameters are a, 6, c, e, 
the matrix S 2 = (Sfj), where 

S = (Sij)= 

/ a b VI -a 2 - b 2 \ 

c VI - c 2 - e 2 e (18) 

V VI - a 2 - c 2 Vc 2 + e 2 - b 2 Va 2 + b 2 - e 2 J 
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is double stochastic. Using the relations V^j = S^j in Eqs. (fT2|) we find values for and cos (5, e.g. 



s\V = , = :, = \/ 1 — a 2 — 6 2 . and slV = , 6 : 
12 Va 2 + 6 2 13 - 23 ^J^2 

rfl1 -& 2 (a 2 + b 2 ) + c 2 (a 2 + 6 2 ) 2 + e 2 (6 2 (a 2 + 1) + a 2 (a 2 - 1)) 

COS u == 

2 a 6 e Vl - a 2 - 6 2 Va 2 + o 2 - e 2 

In the second case when the independent parameters are a, b, d, f one gets 

^ b =, .g> = Vl-« 2 -6 2 , and - ^ + P = ^ 



(19) 



Va 2 + b 2 " ' J Va 2 + b 2 



(20) 



(3) = (a 2 + fr 2 )(6 2 (l - a 2 - b 2 ) - d 2 (a 2 + b 2 )) + / 2 (a 2 + 6 2 (a 2 + b 2 - 1)) 
2a6/Vl - a 2 - b 2 y/a 2 + b 2 - f 2 



and, similarly, in the third case when used parameters are c, d, g, h 



cos 



(3 ) _ Vl-d 2 -fe 2 (3) _ , (3) _ VI -C 2 -d 2 

512 - ^ 2-c 2 -d 2 -h 2 ' s i3-Vc+d+g+h 1, s 23 - V 2 _ c2 _ rf2 _ fe2 (21) 

*< 3 > - (c 2 (l - c 2 ) - d 2 (l - d 2 ) - g 2 (l - g 2 ) + h 2 (l - h 2 ) + (d 2 g 2 - c 2 h 2 )(c 2 + d 2 + g 2 + h 2 - 2)) / (22) 
(2 Vl - c 2 - d 2 ^l-c 2 - g 2 Vl-d 2 - h 2 y/1 - g 2 - h 2 V ' c 2 + d 2 + g 2 + h 2 - l) 

In the first two cases the mixing angles Sjj depend only on three parameters, such that the only constraint on all four 
parameters is given by — 1 < cos 5 < 1 which we write it in the form 

1 - cos 2 5 = sin 2 S > 

By using the form (|19j) for cos 5 the above condition is equivalent to 

sm = 

-(a 2 + b 2 ) 2 (b 2 (b 2 ~2c 2 (a 2 + 6 2 )) + c 4 (a 2 + b 2 ) 2 + e 4 (l - a 2 ) 2 + 2e 2 b 2 (a 2 - 1) + 2e 2 e 2 (a 2 (a 2 - 1) + b 2 (l + a 2 )) 

4 a 2 6 2 e 2 (l - a 2 - 6 2 )(a 2 + b 2 - e 2 ) ~ ° ( ' 23 ' ) 

On the other hand Branco and Lavoura, |33|, and, respectively, Jarlskog and Stora, [32], found a geometric invariant, the area 
of any unitarity triangle, as a measure of CP violation. This area has to be a positive number if the CP symmetry is violated. 
If U, i=l,2,3, denote the lengths of the triangle sides, from (|18[) one gets 



h = ab, l 2 = c\/i - c 2 - e 2 , l 3 = \/l - a 2 - c 2 Vc 2 + e 2 - b 2 (24) 



One makes use of Heron's formula yl = \Jp(j) — h)(p — h)(p — h), where p — (h + h + h)/2 is the semiperimeter, for getting 
area A, and the reality condition A 2 > has the form 

16 A 2 = -(b 2 (b 2 - 2c 2 (a 2 + b 2 )) + c 4 {a 2 + b 2 ) 2 + e 4 (l - a 2 ) 2 +2e 2 b 2 (a 2 - 1) + 2c 2 e 2 {a 2 (a 2 - 1) + b 2 (l + a 2 ))) > 0(25) 



Looking at equations (|23|) and (|25|) we see that they are 
equivalent, the supplementary factors appearing in (|23[) 
being positive for a double stochastic matrix. Hence the 
condition (|16[) on cos (5 is a novel separation criterion be- 
tween the unistochastic matrices and the double stochas- 
tic ones. The separation border is given equivalently by 
A = 0, or cos 2 5=1. 



Equivalent conditions to (|2"51 were given in [31j-[33j. 
where positivity of sin 2 S > was mainly used to con- 
strain the four implied moduli. One aim of the paper is 
to put together all that information for devising a fitting 
model which will use all the possible constraints com- 
ing from unitarity. The idea behind it is to use the non 
uniqueness of the set of four independent moduli used to 
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parameterize a unitary matrix which leads to many re- 
lationships between them and which over constrain the 
error affected experimental data. A consequence will be 
that unitarity implies a very precise tuning between all 
moduli, and we could get information on the physical 
relevant parameters even from a limited experimental in- 
formation. 

From (fT9")) - (l2"Tj) we see that the expressions defining 
and 5 are quite different. The interesting thing is that 
all the possible forms for and cos 6 take the same 
numerical values when the numerical matrix is double 
stochastic. Let us consider the following matrix 



167 


11 


VTT 


200 


20 


200 


8 


</237 


18 


25 


25 


25 


/ 1603 


VMM 


V19253 


v — 


100 


200 



V 40 



0.835 0.55 0.0166 
0.32 0.6158 0.72 
0.4476 0.5642 0.6938 



(26) 



whose entries are within the error bars from (JTJ, and 
the matrix S* 2 = (£?■) is a double stochastic matrix. 
The numbers entering it have been chosen such that the 
measured parameters, V e ±, V e 2, and V T 3, should be 
around the central values from ([TJ. For this matrix we 
get 



„(0 



'23 



110 



V39989 
144 



0.55, 



.(*) 

'13 



200 



0.0166, 



V39989 



0.72, i = 1,2,3 



(27) 



cos 



5419565144 



2066625V211783 



-5.698, i = 1,2,3(28) 



The numerical results show that all the three angles s y - 
are physical, but 5 is not, i.e. the data matrix S, Eg. (1261) . 
does not come from a unitary matrix. Computing the 
triangle area we find 



A = 0.0107? 



(29) 



i.e. an imaginary number which sends the same signal, 
namely, the above data (|26p are not compatible to the 
existence of a unitary matrix. 

Thus the new unitarity constraint — 1 < cos 5 < 1 puts 
non trivial constraints on the experimental data. From 
an experimental point of view the above matrix is per- 
fectly acceptable, from a theoretical point of view it is 
not; and the non-unitarity of the data can be checked by 
the above method. 

In the case when the numerical matrix (JT|) is double 
stochastic, as the above matrix S, Sjj and cos<5 do not 
depend on the chosen four independent moduli used to 
define it. If the modulus of cos 5 is outside the inter- 
val (—1,1) the data are not compatible to the theoreti- 
cal model and the story ends here. If cos S takes values 



within the interval (—1, 1), the data are compatible with 
the theoretical model, and the values for S and Sij nec- 
essary for the recovery of the unitary matrix are easily 
obtained, e.g. from the relations (|19p. and by introducing 
them in Eq.(j4|) we explicitly recover the unitary matrix. 

From almost any numerical experimental data as ([1]) 
we can form a doubly stochastic one by choosing at our 
will four independent parameters, the other five being 
completely determined by equations similar to (|18|) . In 
this way we obtained the equations (fT9"j) - (j2"Tj) for cos<5. 
By using numerical information from matrix (JTJ) for 
a, b, c, d, e, g, h and /, we get 



cos 



cosS { y = -0.9X5. cos.V 1 ,-' = 0.94; 
cosM 2 ) 



-1.1 i, cosSi^ = -0.553, 



r(2) 



cos 



j(3) 

coso_ = 



0.118, cos(5i 2) = 0.737 
0.05i, cos<$( 3) = 0.018«, 
0.0027?:, 



(30) 



where +, c, — denote values of cos S function obtained 
from central values+cr, central values, and, respectively, 
central values— a. The corresponding values for cos d on 
the three rows are different because the numerical val- 
ues for parameters a, b, c, e, respectively, d,f,g,h do 
not come from the same doubly stochastic matrix, e.g. 
e 7^ V 1 — c 2 — d 2 , or numerically 0.7 ^ 0.737. The above 
numerical computations send a clear signal: unitarity of 
the PMNS matrix is a property shared by all its entries 
implying strong correlations between the numerical val- 
ues taken by Vij , and this unitarity feature has to be well 
understood by people doing phenomenology. The border 
of the physical region, A — 0, or cos S = ±1, can be 
used to find maximal intervals for the unmeasured mod- 
uli. If we consider V e i, V e 2, and or V T 3 fixed to the 
central values one gets V^i £ (0.31, 0.45) which represent 
40% from the 3 a interval given by (JlJ , and respectively, 
V^2 € (0.53, 0.64), i.e. about 35% from the 3 a inter- 
val. If one computes now the same quantities by using 
the non unitary double stochastic matrix (|26|) one finds 
Vpi G (0.37, 0.39), and respectively, € (0.57, 0.59), 
i.e. intervals up to 7 times shorter than in the previ- 
ous case. In conclusion the constraints are stronger for 
matrices which are far from unitary ones. The previous 
computations show also how one can modify the numer- 
ical matrix (|26p for obtaining a unitary compatible one. 

Taking into account that more often than not unitar- 
ity constraints are not satisfied by experimental data as 
the numerical results ([3H]) show we have to devise a fit- 
ting model by supplementing the usual x 2- test with a 
separate piece that should implement the fulfillment of 
unitarity constraints, if we want to get reliable values for 
the interesting physical parameters. 
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III. RECONSTRUCTION OF UNITARY 
MATRICES FROM DATA 

We have seen in the previous section that by using 
numerical values from neutrino data, Eq.([T]), one finds 
values for cos<5 outside the physical range [—1,1], and 
generically all numerical values are different. The ex- 
pressions for cos S are provided by the last four relations 
(1121) . and these formulas have to give the same number 
when comparing theory with experiment, by supposing 
the data come from a unitary matrix. Their explicit form 
depends on the four independent moduli used to parame- 
terize the data, see e.g. (Il"9l)- (j2~l"j) . and since there are 57 
such independent groups one gets 165 different expres- 
sions for cos S. 

Hence if data are compatible to the existence of a uni- 
tary matrix the mixing angles s[™' coming from all 57 
groups, and all 165 phases S^ m \ expressed as functions 
of Vodi have to be (approximately) equal, and these are 
the most general necessary conditions for unitarity; they 
can be written as 



< < 1, 



» 



m, n = 1, . . . , 57, (31) 



cos 



<5 W = cos<5 (j) , i,j= 1,...,165 



The above relations are also satisfied by double stochas- 
tic matrices, as the numerical relations from Eqs. ([2"T|) and 
(]28p show, and the condition that separates unitary ma- 
trices from double stochastic ones is given by relation 
([15]). i.e. -1 < cos 6® < 1. 

As a warning, the content of relations (1311) can be sum- 
marized as follows: unitarity implies strong correlations 
between all numerical values V a j. 

Now we present the necessary input that have to be 
taken into account when doing a fit, the experimental 
data are considered moduli of a unitary matrix. It has to 
include as a separate piece the conditions for a complete 
implementation of unitarity. For that we define a test 
function that has to take into account the double stochas- 
ticity property expressed by the conditions (|13|) and the 
fact that in general the numerical values of data are such 
that cos S depends on the choice of the four independent 
moduli and could take values outside the physical range. 
The proposal for the first piece is 



X l = 5>os^)-cos^)) 2 + J2 E 

i<j oc=e,fj,,T \i— 1,2,3 

+ E ( E K- 1 ) > -1<cos<5«<1 (32) 

a=e,/i,r \ i=l,2,3 / 

This proposal was made in for the case of Cabibbo- 
Kobayashy-Maskawa (CKM) quark matrix, and is our 
contribution to the existing methods of reconstructing 
unitary matrices from experimental data, when data are 



compatible with their existence; it expresses the full con- 
tent of unitarity. We stress again that both conditions 
have to be fulfilled: Xi nas to take small values and cos S 
values have to be physically acceptable, since the relation 
Xi = holds true for both double stochastic and unitary 
matrices. When the experimental data are those from 
the matrix (p} the second piece has the form 



xl = 



E E 

— 1,2,3 a.=e,fj.,r 



(33) 



where V a i = \U a i\ are the minimizing parameters, and 
V a i is the numerical matrix that describes the experimen- 
tal data as in ([1]), and a is the errors matrix associated 
to V a %, in our case V e \ — 0.835, a e \ — 0.045, and so on. 
Thus the simplest test function has the form 



2 2.2 
X =Xl+X2 



(34) 



Of course unitarity constraints Q32p could be imple- 
mented by penalty functions, and the concrete form of 
(f3"4"|) could depend on the quality of the fit, which finally 
will choose between the form (|32|) . or that with penalty 
functions. The formula (|34[) provides a least squares 
method which is similar to that used by the BABAR 
Collaboration [27|, the difference being the term Xi that 
implements the theoretical model (QJ. 

The above form can be easily modified if we have other 
experimental information for some functions depending 
on Sij and 5, or on V a i- The primary data in neutrino 
physics are the transition probabilities between two neu- 
trino flavors v a — > t/g, a, ft = e, /i, r at a distance L, 
whose generic form is 35] 



P{v a -> v p ;E, L) = E U^U^U* m U a3 exp(-z A- 

ij 

= EW^) cosA£ + 5>(W*) sinA*. (35) 



i- j 



'■j 



where 

W* p = U aj U; 3 U* m U^ A?- = K 2 - m1)L/2E v (36) 

and a = e, fi, v, i,j = 1,2,3. From [32|, see also [36| . 
we know that the imaginary part has the form 

^(UajU^U^Upi) = J E ea Pi e 3 u ( 37 ) 

with J the Jarlskog invariant. On sees that the parameter 
involving the CP violation appears in transition prob- 
abilities, hence it can be directly determined from the 
experimental data. By using the above relation one gets 



±2 J(sinA? 2 + sin A,, + sin A 



2:i.< 



(38) 
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The term coming from the real part ^(W 7 ^) depends 
explicitly on the flavors a and /?. When a = fi and /3 = e 
one gets 



Similar formulas have to be obtained for all amplitudes 
entering the experimental data. For example the three- 
flavor survival probability P ee in vacuum |5| is given by 



£«(W*)cosA* fc 



(39) 



i?12 COS A 12 + i?i 3 COS A 31 + i?23 cos A 



23 



where 
R12 



R13 = 



R 



23 



u el u e2 ^R(u^u; 2 ) 

~ ( c 12 c 13 s 12( c 23 — s 13 s 23) 
+Cl 2 Ci 3 C23Si2Sl3S23(Ci2 ~ s\ 2 COS 8) 

U el U e3 U^K(U^) 

-C\ 2 c\ 3 Si 3 S 23 (ci 2 S 13 S 23 + C23S12 COSd 

u e2 u e3 u^mu^) 

c\ 3 si 2 si 3 s 23 (ci 2 c 23 cos 5 - S12S13S23) 



(40) 



P ee = 1 Sill 2^13 — COS 



9i 3 sin 2 20i2 sin 2 A- 21 (42) 



l-2(a 2 +& 2 ) + 2(a 2 + b 2 ) 2 



4 a 2 b 2 sin 2 A 2 X 



The above formulas have been obtained by using the 
standard form (j4]) of the PMNS matrix. We use now 
the phenomenological model (|40[) to write Rij in terms 
of invariant quantities, i.e. moduli. One gets 



R12 
R13 



2 2 1 1.2 2 , 2 2 2\ 

a c + be + e — a e ) 



2 2 
a c 



2„2 



b z c 



2 1 2 2 \ 
e + o e J 



(41) 



i? 23 = l(_ 6 2 +a 2 c 2 +5 2 c 2_ e 2 + a 2 e 2 +26 2 e 2 ) 



Hence the transition probabilites, (|40|) -(j42 |) . and 
those similar to them, depends on eight parameters: 
a, b, c, e, J, A12, Ai3,A23 which can be used to ht the 
data. In the usual assumption that U is unitary only six 
are independent, for example a, b, c, J, A\ 2 , A13. That 
means that we may eliminate modulus e. For that we use 
the relation J = 2 A where, J is the Jarlskog invariant, 
and A is the area of (any) unitarity triangle. One gets 



b 2 (l 



c 2 )+a 2 c 2 (l 



{l-a 2 f 



b 2 ) 2e^/a 2 b 2 c 2 (l 



c 2 )(l 



b 2 ) - (1 - a 2 ) 2 J 2 



{l-a 2 f 



r 



(43) 



where e — ±1. If the matrix is unitary then the following 
condition has to be satisfied < e 2 < 1. The second 
parameter is eliminated from the relation 



A? 2 



AI3 



A34 







If we use the the phenomenological models (|2U)) and (f2~T|) , 
i.e. if we take as parameters four independent moduli, we 
find form different terms for all Rij appearing in formulae 



40]) . see Eqs. (gT]) and (|42J). Thus \ 2 could have the form 



EEE 



a P j 



(44) 



where a runs over all transitions probabilities involved 
in the fit, i runs over the sets of four independent mod- 
uli, such as those given by relations (fll?)) - ([2"Tj) . which give 
equivalent but form different expressions for the theoret- 
ical probabilities P{v a ^. Ufl ), and j runs over the number 



of experimental points for each transition probability. By 
doing a fit with the form (I3"4l where \\ is given by (|4"4")) 
one gets values for moduli V a i, Ajj and J. 

The fit could provide a few (approximate) unitary ma- 
trices compatible with the data. In such a case the dou- 
ble stochastic matrices obtained from fit could be used 
to do statistics on unitary matrices. At this point the 
embedding property of unitary matrices into the convex 
set of double stochastic matrices is essential. Indeed if 
U\, U 2 , ... U n are unitary matrices, then 



M 2 = 

i=l 

0<x,< 1, 



i—n 

i=l 
i=l,.. 



1, 



.,n 



(45) 



is double stochastic. The above relation shows that the 
statistics on unitary matrices has to be done through the 
statistics on double stochastic matrices, i.e. the relevant 
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quantities for doing statistics on unitary matrices are the 
squared moduli, |?7aj| 2 . This property allows us to define 
correctly the mean value, (M) , and the error matrix, <jm , 
for a set of doubly stochastic matrices, coming from a set 
of approximate unitary matrices, as follows 



(M) 



\ 



E w /» 



\k=l 



&m = 



\ 



fh—n 



E l^l 4 l n 



(46) 



(47) 



vfc=l 



If the moduli of the mean value matrix (M) obtained in 
this way are compatible with a unitary matrix, or not too 
far from moduli coming from a unitary one, one can re- 
construct from it an (approximate) unitary matrix, and 
in the following we use this technique to reconstruct uni- 
tary matrices compatible to the data, and for obtaining 
information on relevant physical parameters. 



IV. GLOBAL FITS 

By global fit we understand a fit that takes into ac- 
count all available experimental data, and (all) unitarity 



constraints expressed by invariant functions; on this last 
requirement see, e.g., Ref. [37j | . the aim in view being a 
good estimation of mixing and CP violation in neutrino 
physics. For doing that we use the simple ansatz (j3"2")l - 
(|34p in conjunction with data (JlJ to obtain qualitative 
results about J, S and/or V e ^. The last parameter was 
considered free within the interval V e z € (0.02 — 0.4) and 
its error was taken er e 3 = V e g, . Around the central values 
from (JT]) we obtained many approximate double stochas- 
tic matrices, i.e. values for all moduli V a i, and by doing 
statistics with formulae (|47p the following mean value 
matrix and associated error matrix were gotten: 

0.835392 0.537838 0.114933 \ 
Mi = | 0.381844 0.585045 0.715381 (48) 
0.395402 0.607056 0.689253 / 



0.004241 0.00467 0.009417 \ 
<r Ml = | 0.001008 0.00089 0.000263 (49) 
0.003287 0.005701 0.009052 J 

and the unitary matrix corresponding to Mi is 



0.835176 0.537838 0.114933 

Ui= ( -0.0697446- 0.375651 i -0.0445706 + 0.583326 i 0.715381 | (50) 
-0.0668811 + 0.38991 i -0.0434268 - 0.605467 i 0.689253 



If we define two matrices, let us say, by the relation 

M± = JmJ±3^ i (51) 

we get the following prediction for CP violating phase 
and unitarity triangle area 

6 = (89.98ig||)°, A = (127.44 ± 0.01) x 10~ 4 

corresponding to a 3 a level, and so on. 

The results (|38|) and (fl9l strongly depend on the mod- 
uli values around one looks for a minimum. Changing 
these values within the error bars in ([1]) one gets other 
matrices and other values for 6. For example the matrix 

/ 0.851242 0.523442 0.0373568 \ 
M 2 = 0.34007 0.578624 0.741314 (52) 
\ 0.399674 0.625462 0.670118 / 

leads to a value S = (120.01 ± 0.16)°, J 2. x 10~ 4 , 
and so on. This happens since the errors entering ([T]) 
are too large, such that for many choices of fixed three 



moduli, S is practically unconstrained, i.e. there is a 
continuum of unitary matrices with cos S varying in many 
cases between (—1, 1). 

The result (|48|l is "confirmed" by a "standard" fit when 
one uses formulas (|12p , without any approximations com- 
ing from small parameters, etc., and one takes as inde- 
pendent parameters 6ij and 8. The term containing V e 3 
was not introduced in the fit. One gets easily 

0i2 = 0.572332, 6»i3 = 0.088 1 98, (9 23 = 0.79 6 5 47, 
6 = 90.51°, A = 9.9 x 10~ 3 (53) 

where the values for 9ij are given in radians, and \ 2 = 
8.7 x 10~ 2 , where from one gets the following double 
stochastic matrix 

/ 0.837373 0.539489 0.0880842 \ 
M 3 = 0.381844 0.589052 0.712167 (54) 
\ 0.391116 0.601639 0.696462 J 

The corresponding unitary matrix is given by 
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0.837373 0.539489 0.0880842 

U 2 = ( -0.0495463- 0.378656 i -0.0393741 + 0.587735 i 0.712167 | (55) 
-0.0552421 + 0.387195 i -0.0279691 - 0.600988 i 0.696462 



The remarkable fact is that the fit provides a numerical 
value for V e 3 which could not be too far from the true one, 
although the global fit @ failed to find it. To estimate 
the errors to the above parameters we used the moduli 
matrices \JV 2 — 2 a 2 and \JV 2 + 4a 2 , where V and a 
are those given by (Q]), and doing again the fit we got 
X2a = 1-32, and, respectively, xlo- = 1-22, where the 
following prediction was obtained 

8 = (90.51™)°, A = (9.2±§;?) x 10~ 3 (56) 

To see the difference between statistics done by using 
(|47|) and ((STj) and the naive one, we used the central ma- 
trix |T]), and the corresponding error matrix a obtaining 
V c + a, to find 

012 = 0.57996, 0i3 = 0.00054, 

023 = 0.797221, 5 = 79.54° (57) 

results which are equivalent with the moduli matrix 

/ 0.836484 0.547991 0.00054 \ 
W= 0.382939 0.584411 0.715417 (58) 
\ 0.391984 0.598473 0.698697 / 

This matrix is not compatible with a unitary matrix if 
we use our criteria. For example one gets cos^ 3 -* = 
-0.1046 i, and A = 9.0 x 10 -5 i. In fact there are 20 cos 8 
from a total of 165 which take imaginary values, so, in 
principle, the problem cannot be neglected. On the other 
hand one finds x 2 = 6.4 for 4 d.o.f., which is not too bad. 
Taking into account that the double stochasticity prop- 
erty (unitarity) is only weakly violated, we find from (|57j) 
a conservative estimate 8 > 70°. Doing a similar fit for 
V - a one gets \ 2 = 8.1 for 4 d.o.f., and 8 = 91.81°. 

Of course in true fits one does not work with quanti- 
ties such as V ± a, and they were given as pedagogical 
examples. They show that it is possible when using as 
independent parameters 6ij and S to obtain physical val- 
ues for all of them, although the area of unitarity triangle 
could take imaginary values. The above computations 
show that Sjj and 8 are convention dependent, even if 
their numerical values are physical, and in order to min- 
imize this bias we expressed these parameters in terms 
of the moduli V a i which have to be the same irrespec- 
tive of the used convention for the unitary matrix (j4]), 
and we required the fulfillment of all constraints ([3~T]) . 
The second lesson is that one could get values for 8 even 
if moduli where 8 enters are not yet measured, and this 
happens since unitarity over constrains the data implying 
many relationships between the moduli. Of course a neg- 
ative attitude to such results is compliant so long as the 
subtleties implied by unitarity are not fully understood. 



In the following we want to stress the non uniqueness of 
unitary matrices compatible with the error bars in ([1}. 
For example, one can find also a unitary matrix which 
lead to 8 = (75.54 ± 0.11)°; its moduli are 

/ 0.813246 0.575234 0.087959 \ 
Mi = 0.453777 0.614523 0.645328 (59) 
\ 0.364304 0.539877 0.758825 / 

With the matrices Mj, i = 2,3,4, by using the prop- 
erty of double stochastic matrices we find a continuum 
of (approximate) unitary matrices whose squared moduli 
are given by 

M 2 = x Ml + yM 3 + (1 - x - y) M|, < x, y < 1 (60) 

Doing statistics on M,, i = 2,3, 4, with formulae (14"T|) one 
gets 

/ 0.834101 0.546484 0.0750359 \ 
M 5 = 0.394701 0.594258 0.700756 (61) 
\ 0.385327 0.590096 0.694424 J 



( 0.026107 0.023843 0.002995 \ 
om 6 = 0.037521 0.0180207 0.05552 (62) 
\ 0.01148 0.041863 0.06711 / 

With matrix M one can obtain any value for 8 € (75° — 
120°), and, respectively, A € (2 x 10~ 4 - 1.6 x 10" 3 ). The 
prediction at the symmetric point x = y = 1/3 obtained 
from M5 is 

8 = (88.13 twllT, A = (H-1 ± 0-3) x 10~ 4 

at one sigma level given by <7j\/ 5 . 

In the following we want to test the dependence of 
moduli matrix on the form of the unitary matrix U . For 
that we use the Kobayashi-Maskawa form [26| 

(Cl C3S1 S1S3 \ 

C2S1 — C1C2C3 + S2S3 e* v — C1C2S3 - C3S2 e lv J (63) 
sis 2 -C1S2C3 - C2S3 e** -C1S2S3 + C2C3 e 1 ^ J 

where c 2 ; = cos 9i and Sj = sin0i, and ip is the Dirac 
phase. Of course selecting one or another form has no 
theoretical significance because all choices are mathemat- 
ically equivalent; however a clever choice may shed some 
light on important issues. Such issues could be the \U e 3\ 
magnitude and the CP phase 8, (ip), and the above form 
(1531 could be more helpful. 



11 



Doing a fit by using the relations (IT2t and the new 
parameters 0j one gets the moduli matrix (IMf with the 
same digits. In fact we get 

6 l = 0.578337, 9 2 = 0.797342, (64) 
3 = 0.161845, ip = 90.15° 

that confirm the previous results, the only change be- 
ing the value for 6 '3 which is brought, in principle, to a 
measurable value. 

A few conclusions of the above numerical fits could be: 
(a) there are large corridors for moduli around the central 
values in ([]} where data are compatible to the existence 
of a unitary matrix with d ^ 0, tt, and by consequence 
compatible with a CP violation; (b) the real parts of the 
complex entries in unitary matrices are smaller by an or- 
der of magnitude than the imaginary parts, see e.g. ma- 
trix U\, but their neglect completely spoils the unitarity 
property; (c) by consequence, reliable results can be ob- 
tained then and only then when all theoretical formulas 
used in fits are exact, not approximate. 

We are aware that the "experimental" data ([1} are not 
numbers obtained directly from experiments. In the fol- 
lowing we suggest a new kind of global fit able to resolve 



the existence or the lack of the CP nonconservation in the 
leptonic sector by measuring directly the Jarlskog invari- 
ant J, and by testing also the possible new physics which 
could explain the negative value for | U e ^, \ 2 obtained in Q , 
by using the presentday experimental data. For doing 
this we suggest the use of moduli, a, b, c, d, e, /, g, h and 
J, A12, A13, as parameters, although not all of them are 
truly independent. However the moduli are correlated 
through different forms for J and/or cos<5, and relations 
as ([3"T]) . or J(a,h,c,e) = J(c,d,g,h) have to be satisfied 
to a great accuracy. And such a fit could provide val- 
ues for all moduli. If the new fit gives a positive value 
for J one can say that the fit confirms the CP violation 
in the leptonic sector. If the fit will confirm the results 
from Ref. see its Appendix A, concerning the best fit 
obtained for negative values of |t/ e 3| 2 , which imply that 
a and b entering the formula (|42l) for P ee are such that 
1 — a 2 — b 2 < 0, one can use the unitarity constraints 
to obtain other correlations. For example the relation 
(|4*3"f implies also that 1 — a 2 — c 2 = \U T i\ 2 < 0, because 
otherwise e 2 gets complex. If this is true then by using 
as independent parameters a, c, e, and /, and eliminating 
again the parameter e, one gets 



(a 2 + c 2 ){a 2 -c 2 f 2 )-(a 2 -c 2 )f 2 2t^a 2 c 2 f 2 {l ~ a 2 ~ c 2 ){a 2 + c 2 - f 2 ) - (a 2 + c 2 ) 2 J 2 



{a 2 + c 2 ) 2 



(a 2 + c 2 ) 2 



(65) 



Because |J7 r i| 2 = 1 — a 2 — c 2 < 0, the preceding relation 
tell us that It/^l 2 = a 2 + c 2 — / 2 < implying a new 
physics; in fact a new type of unitarity. In any case by 
using the proposed fit there is a bigger flexibility to test 
various consequences of the theoretical model (Q}. 

Of course results such as \U T i\ 2 < 0, i = 1, 2, or 
\U e 3\ 2 < could be artifacts of approximations used for 
the transition probabilities by neglecting the sub leading 
terms, and before speaking of new physics one must use 
exact theoretical forms for all observables measured in 
neutrino physics. The unitarity is a very powerful tool 
because it implies many constraints between the moduli, 
and even a moduli variation of the order of 2-3% could 
spoils it. See in this respect the numerical matrix (f26|) . 
and/or the matrix (|58p. obtained from physical values 
for Sij and 5, which show easily how one can walk into 
a trap. Such a fit will provide values for all moduli even 
if the nowadays data come only from three sources, and 
it deserves to be done since it will give estimates for all 
transition probabilities necessary for the design of future 
neutrino factories. However, if such a result is confirmed 
this will imply a new physics which could be interpreted 
as the first sign showing an evidence for a new generation 
of leptons, or the change of the theoretical frame (j3|, by 
taking the matrix U as being unitary, for example, in a 



non-Euclidean metric. 



V. CONCLUSIONS 

In this work we have discussed the question of the rel- 
evance of the nowadays experimental neutrino data on 
the determination of CP nonconservation in the leptonic 
sector, as well as on solving the negativity problem of 
I C/ e 3 1 2 < 0, raised by the best fit in Ref. @. For doing 
that we used another consequence of unitarity property 
and defined a phenomenological model by taking as free 
parameters the moduli of the unitary matrix, which are 
invariant quantities, and which, naturally, leads to the 
separation criterion between the double stochastic ma- 
trices and those arising from unitary ones. This crite- 
rion provided the necessary and sufficient conditions for 
the consistency of experimental data with the theoreti- 
cal model encoded by the PMNS matrix, the strongest 
condition being — 1 < cos<5 < 1. 

We constructed an exact double stochastic matrix, 
which from an experimental point of view is perfectly ac- 
ceptable, but which does not come from a unitary matrix, 
showing the necessity of using the separation criterion in 
phenomenological analyzes in order to obtain physically 
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meaningful quantities. 

Taking into account the above result, we provided a re- 
liable algorithm for the reconstruction of a unitary ma- 
trix from experimental data, when they are consistent 
with the theoretical model. The fitting algorithm has 
the form a least squares % 2 -test and have two separate 
pieces: one which implements the unitarity constraints, 
and the other which properly takes into account the er- 
ror affected data. We have also shown how this algorithm 
can be modified to take into account the oscillation prob- 
abilities that are the primary data measured in neutrino 
physics. 

We have used the x 2 -test to obtain information on the 
four parameters entering the standard parametrization 
of the PMNS matrix. The test works very well, and by 
applying it we obtained a continuum of (approximate) 
unitary matrices compatible with the nowadays available 
experimental data summarized in the numerical matrix 
([1]). Their explicit construction shows that the moduli 
data ([lj are compatible with the existence of CP viola- 
tion in the leptonic sector, even if one cannot yet find a 
precise value for phase S. 

The numerical results show that unitarity is a very 
strong property, and the tuning of moduli implied by 
it is given at a higher level than expected, the statistical 
errors generated by the approximate character of data are 
at least an order of magnitude less than the experimental 
errors. 



It was also shown that the fitting method, using as 
independent parameters moduli and Jarlskog invariant, 
can directly check the existence of CP violation in the 
presentday neutrino data. On the other hand, the same 
method can resolve the positivity problem for | t/e3 1 2 - If 
the results by Maltoni et al. |5| are confirmed, they will 
be the first signal for a new physics in leptonic sector. 

Last but not least we used the natural embedding (O 
of unitary matrices into the double stochastic set for de- 
vising a method for doing statistics on unitary matrices, 
problem that was an open one until now. The method 
allowed us to find a procedure to obtain mean values, 
and, correspondingly, error matrices starting with a set 
of (approximate) unitary ones. 

As a final conclusion, a satisfactory solution of CP 
nonconservation can be obtained by diminishing the er- 
rors in all running experiments, as well as in all future 
ones, the present theoretical approach being able to test 
the validity of the nowadays unitary model, or the neces- 
sity to modify it. 
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